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ABSTRACT 


The  long  range  Maritime  Patrol  mission  has  evolved  and 
expanded  such  that  the  U.  S.  Navy's  maritime  patrol  aircraft  (P-3C 
Orion)  has  become  an  attractive  platform  from  which  to  employ  a 
wide  variety  of  air  launched  weapons.  Specifically  the  need  for  a 
stand  alone  air-to-air  defensive  capability  was  identified.  In  1989 
the  Naval  Air  Test  Center  at  Patuxent  River,  Maryland  began 
investigation  for  this  through  the  P-3C/AIM-9  integration  program. 
Naval  Postgraduate  School  responded  to  this  with  the  construction 
of  a  ground  vibration  test  stand  and  the  conduction  of  a 
preliminary  vibrational  characterization  of  the  AIM-9  missile.  From 
these  tests  a  two  degree  of  freedom  lumped  mass  model  was  developed 
along  with  the  determination  of  the  primary  and  secondary  missile 
resonant  modes  in  pitch.  In  addition,  a  mathematical  model  of  the 
AIM-9  missile  was  developed  using  finite  element  techniques.  This 
model  was  used  to  analytically  determine  the  modal  parameters  of 
the  missile  and  set  up  a  modal  test  system  to  experimentally  verify 
the  model  parameters  in  view  of  natural  frequencies,  mode  shapes 
and  transient  response.  This  investigation  extends  the  finite 
element  model  utility  by  refining  the  model  structural  code  to 
incorporate  an  improved  mass  distribution  geometry,  thus  providing 
more  accurate  correspondence  between  observed  natural  frequencies 
and  model  generated  modes  and  performing  a  parametric  study  on  the 
model  to  discover  its  sensitivity  to  conformation  and  material 
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I.  INTRODUCTION 
A.  JUSTIFICATION 

As  the  world  political  and  economic  scenes  continue  to  evolve 
and  unfold,  several  facts  seem  to  surface  and  resurface  in  a 
recurring  presentation  throughout  the  theme  of  this  post  cold  war 
environment.  The  first  of  these  is  that  the  United  States  is  the 
only  true  superpower  remaining  on  the  globe.  The  great  majority  of 
the  planet  will  be  looking  to  us  in  the  years  to  come  for  political 
and  economic  leadership,  technological  initiatives  as  well  as  world 
defense  and  regional  security  policy.  Concomitant  with  these 
expectations  is  the  responsibility  of  the  United  States  to  meet 
these  needs.  The  second  fact  is  that  fifty  some  odd  years  of  the 
cold  war  has  left  the  United  States  on  the  brink  of  bankruptcy  in 
an  increasingly  costly  world.  With  the  downfall  of  the  Eastern  Bloc 
nations  along  with  the  Soviet  Union  we  find  ourselves  burdened  with 
supporting  a  much  over-strengthed  and  vastly  miss-apportioned 
defense  structure.  Along  with  this  realization  has  come  an 
increasingly  robust  military  budget  reduction  and  cost  saving 
policy. 

Throughout  our  defense  structure  we  are  witnessing  a  current 
25%  across  the  board  force  reduction  which  promises  to  grow  to  4  0 
or  even  50%  as  Eastern  Europe  and  the  Soviet  Union  become  less  and 
less  potent  militarily.  As  our  nation  presses  on  into  the  1990 's  we 
in  the  Department  of  Defense  are  going  to  be  forced  to  meet  our 


objectives  with  less  material,  a  smaller  force  and  with  less  money. 
The  most  expedient  way  to  accomplish  these  objectives  within  the 
above  constraints  is  to  expand  and  utilize  our  organic  resources 
within  the  research  and  development  purview,  incorporating  into 
these  resources,  to  the  maximum  extent  possible,  all  available 
technological  advances  and  analytical  aids  which  will  develop  the 
most  streamlined  and  efficient  methodologies  possible.  Once  these 
resources  are  built,  incorporated  and  utilized,  hundreds  of 
projects  previously  farmed  out  to  defense  contractors  will  be 
completed  reliably,  in-house,  at  a  savings  of  hundreds  of  millions 
of  dollars. 

In  this  vein  this  investigation  continues  the  ongoing 
development  of  just  such  an  organic  resource. 

B.  BACKGROUND 

Since  the  early  to  middle  1980 's  the  P-3C  Orion  has  employed 
to  great  tactical  effectiveness  the  harpoon  anti-ship  missile 
system.  However,  this  weapon  addition  changed  forever  the  perceived 
threat  potential  to  opponents  as  posed  by  the  P3-C.  Traditionally 
the  Orion  has  been  utilized  and  accounted  as  an  open  ocean  long 
range  anti-submarine  platform.  It's  historical  theater  of  operation 
has  been  mid-ocean  submarine  interdiction,  out  of  the  reach  of  land 
or  sea  based  air  assets.  With  the  development  and  employment  of  the 
sizeable  Soviet  blue  water  Navy,  as  well  as  the  appearance  of  the 
many  third  military  threats  into  the  various  theaters  in  which  the 
P3-C  is  utilized  has  necessitated  the  development  and  incorporation 


of  an  MPA  stand  alone  air-to-air  defensive  capability.  This 
requirement  was  registered  and  quantified  at  the  Naval  Air  Test 
Center  (NATC)  at  Patuxent  River  Naval  Air  Station  under  the 
direction  of  NavAir.  NATC  was  directed  to  develop  and  test  an  air- 
to-air  defensive  combat  system  incorporating  the  AIM-9  Sidewinder 
missile  flown  and  fired  from  the  P3-0RI0N.  This  developmental 
project  frontier  a  quantum  leap  in  mission,  tactics  and 
capabilities  within  the  Maritime  Patrol  community.  For  the  first 
time  the  Orion  would  have  the  capability  to  encounter,  and  with  a 
reasonable  chance  of  success,  defeat  an  air  threat.  The  Harpoon 
anti-ship  mission  combined  with  this  new  anti-air  capability  placed 
the  P3-C  high  on  the  potential  enemies  air  threat  priority. 

There  are  many  factors  which  must  be  considered  when 
attempting  to  employ  a  new  weapon  system  from  an  established 
platform.  Primary  amongst  these  within  the  air  community  is  the 
static  and  dynamic  interaction  between  the  wing  and  the  proposed 
weapon  store.  Weapon  induced  wing  flutter  is  the  first  and  most 
lethal  aerodynamic  phenomenon  to  be  addressed. 

Every  structure  or  entity  that  possesses  mass  and  occupies 
space  has  as  physical  characteristics,  natural  modes  and 
frequencies  at  which  they  will  vibrate  under  excitation.  When  a 
structure  (or  conglomerate  of  structures)  is  excited  at  its  natural 
frequency,  great  increases  in  vibrational  amplitude,  leading  to 
structural  failure,  can  occur.  Such  was  the  case  for  the  Tacoma 
Narrows  bridge.  When  excited  by  a  forcing  function  (in  this  case 
the  wind)   at  the  spans  natural   frequency,   the  vibrational 


amplitudes  gradually  increased  until  catastrophic  failure  of  the 
bridge  occurred.  This  case  was  uniquely  noteworthy  because  the 
increases  in  vibrational  amplitude  occurred  over  the  course  of 
several  hours  and  was  well  documented.  However,  the  catastrophic 
failure  of  an  aircraft  wing  due  to  flutter  induced  by  exceeding 
design  performance  specifications  or  wing  to  wing  store  interaction 
usually  occurs  within  very  few  cycles  and  rarely  exceeds  more  than 
one  or  two  seconds  to  failure  from  onset.  Thus  the  utility  and 
desirability  of  an  accurate  flutter  prediction  methodology  is  self 
manifest. 

C.  PREVIOUS  WORK 

This  thesis  is  the  third  in  an  ongoing  investigation  to 
design,  build  and  implement  a  Navy  organic  ground  vibration  test 
facility  and  wing  flutter  prediction  methodology. 

Initially  Lieutenant  Commander  J.  B.  Hollyer  designed  and 
built  a  Ground  Vibration  Test  stand  (GVT)  and  utilized  it  to 
ascertain  experimentally  the  first  three  modes  in  pitch  of  the  AIM- 
9  Sidewinder  missile.  From  this  work  a  two  degree  of  freedom  lumped 
mass  model  was  developed  to  represent  the  missile  dynamics. 
Utilizing  Lagrange's  equations  the  following  is  presented  in  review 
for  the  lumped  mass  model.  Figure  1  presents  the  spatial  geometry, 
the  coordinate  systems  used  and  the  component  make-up  of  the  lumped 
mass  model.  Figure  2  depicts  the  various  coordinate  assignments 
used.  This  is  followed  by  the  assumptions  made  for  the  theory  and 
organizational  definitions  for  both  the  mass  energy  equations  and 
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Mass  Spring  Damper  System 
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Coordinate  System 
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the  State  Space  format. 
Assumptions; 


sin  e=e 


Mass  connections  are  rigid: 

M!  =  M2  =  M3  =  M 

Definitions: 


(1) 


Cos  9=1  (2) 


yx=y+l*&  o) 

Y^Y  (4) 

Y2=Y-L*Q  (5) 


T,   =  Kinetic  Energy,   =  .  ,  dY,  v  2  ,-. 

1  yjri  1/2  *#.*(—!)  (6) 

1   dt 


d7  2 
T,  =  l/2*Af?*(— -^)  (7) 

2       2    dt 


dY     2 
T3  =  l/2*Ms*(—Z.)  (8) 

J      J   dt 


V  =  Potential  Energy  = 


(9) 

Mt*g*  Y±+M2* g*  Y2+M3*g*  Y^+l / 2  *KX*  (x+L/ 2  *6)  2  +1/2  *K2*  (x- 1/2  *0)  2 


Using  Laaranaian  Equations; 

L'=T-V  (10) 


d   .     dL'    v  _  dL1  _Q 

dt     ^dq,        dq±    Ui  (11) 

o — —- 
dt 


T=l/2*(M1*(^+L*d^)2+M2*(^)2+M3*{-d£-L*d-%-)2)  (12) 

dt  dt         z      dt         *      dt  dt 

q,  @  Y: 

J^=l/2*(M*(2*(-^+L*dA)+2*M*-^+2*M*(-^-L*dA)))=M*(-^+L* 
oag/!  at     dt  dt  dt  dt  dt 

dt  dt  dt  dt 

(13) 


J*         *jf       =3*M*d2^- 

dt  a(dq)  dt2  (I4) 

dt 
q2  @  6: 


=M*L*(^+L*-^)+M*L*(-^-L*-dQ)=(2*M)  *L*dY 


_.dg2N       dt         dt  dt         dt  dt       (15) 

a("dF) 


a(-dF' 

V  =  Potential  Energy= 

M*g*(Y+L*Q)  +M*g*Y+M*g*(YL*Q+l/2*K1*{Y+L/2*B)2+l/2*K2*(Y-L/2*e)2 

(17) 


q,   @   Yi 


7  =3*M*g+K*(Y+  —  *9)  +tf*(r— £*©)  (18) 


3g1  2  2 


q2   @   0: 


'=iq*^*(y+^*e)-iC2*-£*(y-i:*e)  (19) 


dg2      1     2  2  2     2  2 


Q.*: 


3*M*-^-^+(M1-Mj  *L*-^-5+2*J:*y+(Ji:,-Jl)  *Le+3*M*sr=F,  (20) 

dt2         1      3  dt2 


Q 


2#' 


2*M*L*-^-^+(M1-Mj  *L2*-^-^  +  (X1+/C)  *  —  *r+(JC-^)  *—  *Q+(Mr-M%)  *g*L=Fl 

dt2  dt2  2  x      2         4  *      3  * 

(21) 


# 


3*M*^-^+1'3(K1+K2)  *7=F,  (22) 

dt2 


2*M*L*^-X+2'*2*K*—  *Y  (23) 

dt2         2 


l=(C1  +  C2)*f 


(24) 


2=<C1  +  C2)*|*f 

4=(C1-C2)*^*d-| 

Dropping  the  last  term  on  the  left  sides  of  Equations  23  and  22  due 
to  unnecessary  initial  conditions  and  proceeding  into  the  state 
space  format  yields: 

Y=X,       e=x3 
dY/dt=X2       de/dt=X4 
dX,/dt=  X2 

dX3/dt=  X4 


X  =    [Xx  X2  X3  XA]T  (25) 

dX 
3*M*—±=-2*K*X1-(Cl+C2)  *X2-(C1-C2)  *X<  (26) 


2*M*L*— -2-=-2*JO<-*X1-(C1  +  C2)  *—*X9-(C,-C9)  *—*X.  (27) 

dt  2       x         i2221244 


In  = 


10  0      0 

0  (M1=M2+M3)  0  (M1-M2)  *L 

0      0  10 

0  (M1+M2)*L  0  (Afj-ilfj)  *L2 


(28) 


An= 


0  1 

~{kx-k2)  -(q+c2) 

0  0 

(*1+ie2)*!  -(q+c2) 


<*i-*2)*-f 


(jq-jc,)*- 


-(q-c2) 

i 


(29) 


The  homogeneous  problem,  representing  a  system  response  to  an 
initial  condition  is: 


n    dt      n 


(30) 


Premultiply  by  [!„]"'  to  obtain: 


dt     n        n 


The  eigenvalue  problem  can  be  formed  as 


(31) 


[  X*I-A  ]  *X=[Q] 


(32) 


For  a  non-trivial  solution,   one  obtains  the  charachteristic 
equation 


\X.*I-A\  =  0 


(33) 
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Once  the  eigenvalues,  X;  (i=l  to  n)  has  been  found,  the  mode  shapes 
may  be  identified  using  the  relation  that: 


\.*XU)    =  A*XU)  (34) 


where  X(i)is  the  state  vector  corresponding  to  the  i*   eigenvalue  \t. 

X(t)=eMt*X(0)  (35) 

provides  the  solution  for 

S!1£=A*X+B*UU)  (36) 

After  the  construction  of  the  GVT  and  the  determination  of  the 
principal  modes  in  pitch  by  Lcdr.  Hollyer,  finite  element 
techniques  were  used  by  Mr.  Michael  Shutty  to  construct  an  initial 
finite  element  model  of  the  AIM-9  Sidewinder.  Utilizing  this 
prototype  model  Mr.  Shutty  developed  a  methodology  for  utilizing 
and  exploring  the  various  applications  of  the  PAL-2  software  used 
to  construct  the  finite  element  model. 

D.  PURPOSE 

The  purpose  of  this  thesis  was  to  further  the  development  of 
an  existing  but  incomplete  methodology  for  a  Navy  organic  ground 
vibrational  test  and  prediction  capability.  In  order  to  accomplish 
this,  integration  of  prior  work  by  Lcdr  Hollyer  and  Mr.  Shutty  was 
utilized  within  the  framework  of  this  investigation. 
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E.  SCOPE 

The  evaluation  includes  refinement  and  reintegration  of  the 
AIM-9  finite  element  model.  This  incorporates  novel  mass  location 
and  distribution  logic  as  well  as  real  material  characteristics 
reflecting  accurate  representation  of  missile  construction.  This 
logic  was  honed  to  most  accurately  represent  the  missile  in  terms 
of  modal  shapes  and  freguencies  as  determined  experimentally.  A 
parameterization  of  missile  model  properties  was  then  conducted  to 
determine  the  models  sensitivity  to  geometric  as  well  as  material 
changes . 

The  parameterization  primarily  focused  on  missile  component 
mass  and  modulus  of  elasticity  changes  as  well  as  the  replacement 
of  selected  structures  with  potential  future  improvements  in  terms 
of  material  characteristics. 
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H.  INVESTIGATION 

A.  GENERAL 

Finite  element  analysis  (FEA)  is  the  prediction  of  structural 
modes  and  frequencies  utilizing  a  finite  element  model  to  undergo 
computer  generated  excitation.  A  finite  element  model  is  the 
representation  of  a  continuous  structure  by  a  discrete  one.  In  the 
finite  element  method  structures  are  descretized  in  that  they  are 
represented  by  discrete  grid  or  node  points  connected  by  structural 
elements.  This  representation  is  in  the  form  of  a  mathematical 
model  consisting  of  discrete  elements  connecting  discrete  node 
points.  The  rule  of  thumb  applied  to  discretization  and  thus  to  FEA 
is  that  the  more  node  points  used  the  more  accurate  the  solution 
that  is  obtained.  Also  known  is  that  the  fewer  the  nodes  the 
stiff er  the  model.  Unfortunately,  for  every  node  point  incorporated 
into  the  model  six  degrees  of  freedom  are  introduced  into  the 
computational  analysis.  This  rapidly  becomes  the  limiting  factor  in 
terms  of  computing  time  and  space  as  the  coefficient  matrices 
become  unwieldy. 

B.  FINITE  ELEMENT  MODEL 

The  refined  MSC/PAL2  software  code  developed  for  this 
investigation  of  the  AIM-9  Sidewinder  geometry  is  presented  in 
Appendix  A.  As  the  code  illustrates,  729  node  points  were  used 
which  generated  4374  dynamic  degrees  of  freedom  for  the  excited 
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structure.  Figure  3  depicts  missile  geometry  as  generated  by  the 
code  upon  building  the  model. 


Refined  Model  Geometry 


Figure  3 
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The  large  number  of  nodes  was  required  to  mathematically  represent 
the  missile.  Instead  of  using  curved  beam  segments  which  the 
software  does  not  support,  concentric  rings  consisting  of  16  nodes 
joined  with  small  quadrilateral  plate  elements  were  utilized  to 
represent  the  missile  cylindrical  body.  The  concentric  rings  were 
spaced  closely  at  three  inches  apart  to  avoid  aspect  ratios  within 
the  plate  elements  of  greater  than  four  to  one,  as  the  dynamic 
response  becomes  unreliable  above  this  ratio.  The  resulting 
geometry  yielded  thirty  five  concentric  rings  of  node  points 
comprising  the  uniform  cross  section  of  the  cylindrical  body  and 
three  decreasing  radius  rings  of  node  points  comprising  the 
tapering  nose  cone.  Quadrilateral  plate  elements  of  .25  inch 
thickness  were  used  in  defining  the  skin  of  the  cylinder  and  nose 
cone,  while  triangular  plates  of  .25  inch  thickness  were  used  to 
describe  the  cone  tip  and  the  rear  plate  of  the  missile.  The  AIM-9 
attaches  to  the  LAU-7  launcher  rail  at  three  points.  This  was 
modeled  by  zeroing  the  displacements  at  points  85,  213  and  373  in 
the  x,y  and  z  directions  as  well  as  zeroing  the  rotations  about  the 
y  and  z  axes.  Rotation  about  the  x  axis  was  not  zeroed  since  the 
motion  in  pitch  is  of  primary  interest  here. 

The  front  fins  were  constructed  using  a  combination  of 
quadrilateral  and  triangular  plates  each  of  0.10  inch  thickness. 
The  AIM-9  front  fins  are  actually  moveable  canards  which  rotate 
symmetrically  about  their  attachment  axis.  Initially  it  was  thought 
that  this  feature  would  pose  an  attachment  problem  to  accurately 
represent  the  canards  for  dynamic  analysis.  This  however  did  not 
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impact  modeling.  The  concern  of  this  investigation  is  primarily 
developing  a  methodology  for  predicting  wing  flutter  induced  by 
weapon/wing  interaction.  In  order  for  flutter  to  take  place  the 
weapon  must  be  attached  to  the  wing  and  the  missile  maneuvering 
canards  would  be  disabled.  The  canards  were  attached  to  the  missile 
via  a  .035  inch  plate.  Due  to  the  small  size  of  the  attachment 
plate  the  front  fins  exhibited  an  insufficient  stiffness  problem 
which  could  be  addressed  easily  with  further  work  on  a  subsequent 
model . 

The  rear  wings  were  modeled  using  a  combination  of  0.5  inch 
quadrilateral  and  triangular  plates.  These  were  rigidly  attached  to 
the  missile  body  using  the  laterally  arrayed  node  points  of 
successive  concentric  body  rings. 

The  mass  distribution  node  points  represent  the  most 
significant  refinement  of  the  finite  element  model.  Twenty-two  mass 
distribution  nodes  were  located  along  the  geometric  centerline  of 
the  missile.  Each  missile  component  (i.e.  motor,  seeker,  warhead, 
etc.)  received  it's  allotted  number  of  nodes  depending  on  its 
centerline  axis  dimension.  These  nodes  were  then  secured  spatially 
to  the  missile  via  a  wagon  wheel  spoke  arrangement  which  connected 
each  mass  node  to  the  nodes  located  on  the  corresponding  concentric 
cylindrical  body  ring.  These  connection  elements  were  then  assigned 
material  properties  according  to  what  structure  they  were  to 
materially  represent  (i.e.  Propellent,  explosive,  seeker  etc.).  The 
mass  for  each  missile  component  was  then  evenly  distributed  among 
the  appropriate  nodes.  This  arrangement  provided  the  most  accurate 
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representation   of   the   actual   missile   structure   and   weight 
distribution . 

C.  MODE  AND  FREQUENCY  CORRESPONDENCE 

Completion  of  the  model  geometry  enabled  the  pursuit  of  the 
next  investigative  goal,  correspondence  of  the  model  generated 
modes  and  freguencies  to  those  determined  experimentally.  In  order 
to  accomplish  this,  seguential  static  and  dynamic  response 
evaluations  utilizing  the  "Build  Model",  "Statics"  and  "Dynamics" 
features  of  the  software  were  reguired  to  assess  the  affect  of  the 
conformal  adjustments  made  on  the  model.  These  adjustments  were 
performed  on  a  trial  and  error  basis  in  order  to  achieve  the 
highest  degree  of  correspondence  to  the  experimentally  derived 
parameters.  In  using  the  "Build  Model"  feature  the  system 
stiffness,  damping  and  mass  matrices  were  formed.  These  were 
compiled  based  on  Newton's  Second  Law  of  Motion  which  provides  the 
foundation  upon  which  all  vibrational  work  is  performed.  Written  as 
a  second  order  differential  eguation  in  the  time  domain  we  obtain 
the  following: 


Where 


M*Z(t)  +C*Z(t)  +K*Z{t)  =f{c)  (37) 


M  =  system  mass  matrix 

C  =  system  damping  matrix 

K  =  system  stiffness  matrix 

x{t)=  nodal  displacement  vector  (38) 

xU)   =  nodal  velocity  vector 

x{c)   =  nodal  acceleration  vector 


f{t)   =  applied  force  vector 
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The  system  equations  formed  during  a  modes  analysis  include 
the  system  stiffness  and  mass  matrices,  which  when  placed  in  matrix 
notation  yields  the  following: 


M*x   +  K*x   =  0  (39) 

This  function  is  set  to  zero  since  there  is  no  forcing  function  in 
normal  modes  analysis.  All  motion  is  assumed  to  be  steady  state 
harmonic  motion.  The  missile  oscillates  at  some  frequency  f(t)/  with 
each  point  moving  either  in  or  out  of  phase.  Thus  the  node 
acceleration  vector  may  be  written  in  terms  of  the  displacement 
vector  as: 

x   =  -(2*n*f)2*x  (40) 

Substitution  then  yields: 

[  K  -    (2*n*f)2*M  ]  *x  (*D 


Excluding  the  trivial  solution,  the  remaining  non-trivial  solution 
to  equation  (39)  produces  the  system  of  natural  frequencies  for  the 
specific  finite  element  model.  Associated  with  each  of  these  is  a 
companion  mode  shape.  As  there  is  no  forcing  function  utilized 
during  the  normal  modes  analysis,  the  scaling  parameters  are 
independently  arbitrary  so  that  only  the  relative  mode  shape  holds 
any  meaning. 

18 


The  normal  modes  analysis  in  this  study  were  accomplished 
utilizing  the  sub-space  iteration  feature  on  the  software  since 
only  the  lowest  few  resonant  frequencies  were  sought.  These  were 
selected  because  in  reality  only  the  first  few  natural  frequencies 
are  excited  to  resonance  by  the  dynamic  forces  applied  to  the 
system.  Therefore  these  are  of  primary  interest  in  predicting  wing 
flutter.  The  sub-space  iteration  method  solves  for  the  eigenvectors 
and  the  eigenvalues  utilizing  the  mass  and  stiffness  matrices  in 
equation  (40) .  These  matrices  are  derived  from  the  material 
characteristics  coded  into  the  model.  The  results  of  the  static  and 
dynamic  analysis  using  sub-space  iteration  on  the  finite  element 
model  are  presented  in  Appendix  A. 

Evaluation  of  the  refined  models  initial  modal  analysis 
revealed  that  the  three  lowest  fundamental  frequencies  in  pitch, 
which  were  the  frequencies  of  interest  for  this  study,  were 
significantly  different  than  the  experimentally  derived  values. 
This  indicated  that  the  physical  characteristic  parameters  within 
the  model  code  required  adjustment  to  more  accurately  represent  the 
real  missile.  This  was  accomplished  in  increments.  First,  the 
modulus  of  the  propellent  was  adjusted  according  to  published  data 
on  propellent  formulation  and  structuring  criteria.  This  brought 
the  frequencies  into  correct  perspective  relative  to  each  other. 
The  modulus  of  the  warhead  was  changed  next  on  a  trial  and  error 
basis  in  an  attempt  to  bring  the  frequencies  into  closer 
correspondence  with  the  experimentally  derived  values.  This  effort 
was   successful   as   evidenced   by   figure   4.   This   graphical 
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representation  of  the  first  three  natural  frequencies  as  each 
improvement  to  the  finite  element  model  was  made,  depicts  the 
progressive  honing  of  the  model  towards  the  physical  reality  of  the 
actual  missile.  The  experimentally  obtained  values  are  depicted  by 
the  white  diamonds,  and  represent  the  first  three  natural  modes  of 
the  AIM-9  Sidewinder.  The  first  three  modes  as  produced  by  the 
original  finite  element  model  are  shown  by  the  black  squares  with 
grey  circle  centers.  As  is  evident  by  the  graph,  there  is  quite  a 
disparity  in  the  values,  particularly  in  the  second  mode.  With  the 
addition  of  the  center line  axis  for  mass  distribution,  as  shown  by 
the  white  squares  with  grey  centers,  we  see  a  marked  improvement  in 
second  mode  correspondence  with  a  degradation  in  first  and  third 
correspondence.  As  the  propellent  modulus  was  adjusted  an 
improvement  in  first  and  second  modes  is  noted  with  further 
degradation  in  the  third  mode.  Here  for  the  first  time  correct 
relative  position  between  modes  is  noted.  With  the  adjustment  of 
the  warhead  modulus  as  represented  by  the  grey  diamonds,  marked 
improvement  in  all  three  modes  is  seen  for  the  closest 
correspondence  to  the  experimental  values.  Further  improvement  in 
frequency  correspondence  could  be  obtained  with  continued 
refinement  of  the  missile  component  physical  characteristics.  At 
the  conclusion  of  this  investigation  model  frequencies  were  within 
6%  for  modes  1  and  3,  and  within  18%  for  mode  2.  Table  one 
illustrates  the  frequencies  presented  in  Figure  4  as  obtained  in 
the  progressive  refinement  of  the  finite  element  model. 
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Table  1. 


MODE  1 

MODE  2 

MODE  3 

ORIGINAL  MODEL 

29.7  HZ. 

50.8  Hz 

53.0  Hz 

MEASURED  DATA 

24.4  Hz 

34.0  Hz 

52.1  Hz 

1st  CENTERLINE 
RUN 

34.4  HZ 

43.5  Hz 

49.1  Hz 

ADJ.  PROPEL. 
MODULUS 

21.0  Hz 

23.6  Hz 

42.7  HZ 

ADJ.  WARHEAD 
MODULUS 

23.7  HZ 

28.1  Hz 

52.1  Hz 

C.  PARAMETRIC  STUDY 

Versatility  and  adaptability  earmark  a  truly  utile 
methodology.  It  was  identified  as  desirable  while  developing  this 
capability,  to  maintain  to  the  greatest  extent  possible  the  ability 
to  incorporate  and  accommodate  physical  block  changes  in  the 
subject  missile.  Research  and  development  will  yield  more  efficient 
and  durable  propellants,  warheads,  guidance  packages  and  seeker 
heads  to  be  incorporated  into  future  models  of  existing  weapons  or 
totally  new  kill  stores.  These  changes  must  be  accounted  for  to 
accurately  predict  wing-kill  store  interactive  wing  flutter. 
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To  investigate  this  methodology's  sensitivity  to  any  physical 
block  change  in  the  subject  weapon,  a  parametric  study  was 
conducted  using  the  refined  finite  element  model  for  the  AIM-9 
Sidewinder.  Gross  mass  and  modulus  of  elasticity  changes  were  made 
for  the  warhead  and  rocket  motor.  These  two  missile  components  were 
chosen  because  they  comprise  the  majority  (upwards  of  75%)  of  the 
missiles  mass  and  length.  With  these  changes,  evaluations  were  made 
to  determine  the  impact  they  had  on  1.  The  modal  characteristics  of 
the  missile,  and  2.  The  systems  calculation  completion  sensitivity 
(will  the  program  crash) .  Finally,  the  missile  casing  material 
parameters  were  changed  to  represent  a  potential  future  composite 
casing.  Again  an  evaluation  was  performed  to  determine  the  modal 
ramifications  given  composite  casing  replacement. 

In  varying  the  modulus  of  elasticity  of  the  rocket  motor 
propellent  it  was  desired  to  examine  the  ramifications  on  mode 
freguency.  In  anticipation  of  future  improvements  in  propellent 
yield,  burn  temperature  and  density  only  lower  than  present  modulus 
examples  were  evaluated.  Evaluations  at  50%  original  and  25% 
original  modulus  were  performed.  The  data  results  are  presented  in 
Appendix  A. 

The  variation  of  the  warhead  weight  was  accomplished  over  the 
entire  lineal  length  of  the  missile.  As  the  anticipated  weight 
saving  advancements  mentioned  above  are  accomplished,  more 
explosive  may  be  added  to  the  weapon  to  increase  its  effective  kill 
radius.  Thus  the  weight  variation  study  took  place  over  a  range  of 
125%  mass  increase  to  60%  mass  reduction  in  the  warhead  component. 
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The  data  results  are  presented  in  Appendix  A.  Simulation  of  a 
composite  missile  skin  was  desired  in  anticipation  of  its 
incorporation  into  the  actual  production  weapon  in  the  near  future. 
This  would  engender  both  a  significant  weight  savings  on  the  order 
of  20%,  and  provide  a  strength  and  flexibility  improvement  which 
would  translate  directly  into  a  missile  performance  enhancement.  A 
currently  available  carbon  fiber  composite  was  utilized  with  the 
physical  material  data  incorporated  into  the  model  by  changing  the 
material  component  parameters.  The  data  results  are  presented  in 
Appendix  A.  This  simulation  represents  all  of  the  quadrilateral 
plates  comprising  the  missile  skin  as  similar  plates  made  of  the 
chosen  carbon  fiber  composite.  All  of  the  connection  methods  and 
geometric  specifications  remain  the  same. 
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III.  CONCLUSIONS 

In  formulating  this  study  it  quickly  became  apparent  that  the 
tools  necessary  to  thoroughly  accomplish  the  development  of  an 
accurate  wing  to  wing  store  interactive  flutter  prediction 
methodology  were  already  in  existence.  The  wing  stores  to  be 
examined  readily  lent  themselves  structurally  to  finite  element 
modeling  and  an  increasingly  accurate  representation  could  be 
accomplished.  The  key  to  the  more  accurate  solutions  lay  in 
understanding  the  software's  capability  to  geometrically  compose 
the  various  weapon  components  and  the  interrelationships  these 
different  components  made  in  the  formulation  of  the  systems  of 
equations  which  represented  the  physical  dynamic  response.  As  an 
example,  the  incorporation  of  the  centerline  mass  distribution  axis 
provided  an  increased  capability  to  distribute  the  mass  components 
of  the  model  in  a  more  refined  manner.  Thus  an  increased 
correspondence  to  the  experimental  data  in  the  first  three  natural 
frequencies  was  accomplished.  In  the  parametric  study  mass 
increases  or  decreases  didn't  affect  the  modes  as  much  as  where  the 
mass  was  located  spatially.  This  was  evident  in  the  differences 
seen  between  the  original  model,  which  had  the  masses  placed  on  the 
skin  of  the  model,  and  the  refined  centerline  mass  model  which 
distributes  the  mass  evenly  along  the  geometric  centerline. 
Accordingly,  individual  component  modulus  changes  did  not  impact 
modes  as  did  the  relative  positioning  of  different  modulus  values 
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along  the  missile  length.  A  combination  of  mass,  modulus  and 
location  /juxtaposition  of  these  values  could  significantly  impact 
model  natural  frequencies. 

The  techniques  developed  and  described  in  this  investigation 
are  worthy  of  further  work.  The  methodology  accumulating  throughout 
the  three  thesis'  to  date  will  ultimately  result  in  the  ability  to 
accurately  predict  potential  weapon/wing  store  interactive  wing 
flutter  for  the  price  of  a  technicians  time.  This  would  save 
millions  of  taxpayer  dollars. 
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IV.  RECOMMENDATIONS 

Further  refinement  of  the  finite  element  model  is  possible  and 
required.  Continued  refinement  of  the  fin  attachment  points  to  the 
missile  body  would  further  increase  the  realism  of  the  mathematical 
representation.  Stiffening  of  the  existing  attachments  plus  the 
addition  of  multiple  adjacent  attachments  are  potential  solutions. 
The  aft  end  plate  also  requires  some  modification.  The  actual 
missile  incorporates  a  Con-Di  nozzle  especially  engineered  to 
maintain  a  constant  burn  rate  over  the  firing  life  of  the 
propellent.  This  must  be  simulated  in  replacement  of  the  flat 
circular  end  plate  now  incorporated.  This  could  be  accomplished 
utilizing  converging  triangular  plates  into  the  aft  end  of  the 
missile. 

The  GVT  stand  analysis  hardware  could  also  be  upgraded  and 
optimized.  Replacing  the  current  data  accumulation  and  analysis 
equipment  with  a  spectrum  analyzer  would  enable  the  rapid 
experimental  determination  of  the  systems  natural  frequencies 
referenced  to  all  axes  of  interest  with  a  minimum  of  re-program  and 
alternate  hardware  placement. 

Further  investigation  of  the  GVT,  launcher,  rail  interaction 
in  all  axes  of  interest  would  be  required  to  determine  the  amount 
of  interference  (if  any)  being  produced  between  the  weapon  and  the 
wing.  Stiffening  of  the  GVT  would  also  be  required  if  the  two  other 
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primary  axes  of  motion  were  to  be  investigated  as  yaw  and  roll  axes 
signals  are  overcome  with  interference. 
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APPENDIX  A 


MSC/pal  2 


Mode   1 


Figure   5. 
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nSOpal  2 


Mode   2 


Figure   6 
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MSOpal  2 


Mode   3 


Figure   7 
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MODE  and  FREQUENCY  DATA 


ORIGINAL  WEIGHT  DISTRIBUTION  ON  CENTERLINE  AXIS 
12-09-91   14:21  MSC/PAL  2   V  4.0 

-  AIM9  GEOMETRY  FILE 


Page 


MODE  NO.  1  AT  3.43604E+01  CPS  (  2.15893E+02  RAD/SEC  ) 
MODE  NO.  2  AT  4.35337E+01  CPS  (  2.73530E+02  RAD/SEC  ) 
MODE  NO.     3  AT   4.91267E+02  CPS   (  3.08672E+02  RAD/SEC   ) 


MODULUS  AT  .25  ORIGINAL  VALUE 

12-11-91   11:23  MSC/PAL  2   V  4.0 

-  AIM9  GEOMETRY  FILE 


Page 


MODE  NO.  1  AT  2.10428E+01  CPS  (  1.32216E+02  RAD/SEC  ) 
MODE  NO.  2  AT  2.32971E+01  CPS  (  1.46380E+02  RAD/SEC  ) 
MODE  NO.     3  AT  4.24137E+01  CPS   (   2.66493E+02  RAD/SEC   ) 


MODULUS  AT  .50  ORIGINAL  VALUE 

12-10-91   17:53  MSC/PAL  2   V  4.0 

-AIM9  GEOMETRY  FILE 


Page 


MODE  NO.     1  AT  2.10642E+01  CPS   (   1.32350e+02  RAD/SEC   ) 
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MODE  NO.     2  AT  2.33299E+01  CPS   (   1.46586E+02   RAD/SEC  ) 
MODE  NO.     3  AT   4.26267E+01  CPS  (  2.67831E+02  RAS/SEC  ) 


WARHEAD  MASS  AT  .50  ORIGINAL  VALUE 

12-12-91    11:11  MSC/PAL  2  V  4.0         Page 

-AIM9  GEOMETRY  FILE 

MODE  NO.  1  AT  2.10572E+01  CPS  (  1.32306E+02  RAD/SEC  ) 
MODE  NO.  2  AT  2.33381E+01  CPS  (  1.46638E+02  RAD/SEC  ) 
MODE  NO.     3  AT  4.26968E+01  CPS  (   2.698272E+02  RAD/SEC  ) 


WARHEAD  MASS  AT  .40  ORIGINAL  VALUE 

12-04-91    15:17  MSC/PAL  2  V  4.0  Page 

-AIM9  GEOMETRY  FILE 

MODE  NO.      1  AT  3.00438E+01  CPS  (   1.88771E+02  RAD/SEC   ) 

MODE  NO.      2  AT  5.09137E+01  CPS  (   3.19900E+02  RAD/SEC   ) 

MODE  NO.      3  AT  5.29758E+01  CPS  (   3.32857E+02  RAD/SEC   ) 


WARHEAD  MASS  AT  .60  ORIGINAL  VALUE 

12-04-91    17:15  MSC/PAL  2  V.  4.0  Page  2 

-AIM9  GEOMETRY  FILE 
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MODE  NO.  1  AT  3.02458E+01  CPS  (  1.90040E+02  RAD/SEC  ) 
MODE  NO.  2  AT  5.09760E+01  CPS  (  3.20291E+02  RAD/SEC  ) 
MODE  NO.     3  AT  5.29784E+01  CPS  (   3.32873E+02  RAD/SEC   ) 


WARHEAD  MASS  1.2  5  ORIGINAL  VALUE 

12-11-91     18:34  MSC/PAL  2  V  4.0 

-AIM9  GEOMETRY  FILE 


Page  3 


MODE  NO.  1  AT  2.10205E+01  CPS  (  1.32076E+02  RAD/SEC  ) 
MODE  NO.  2  AT  2.32820E+01  CPS  (  1.46285E+02  RAD/SEC  ) 
MODE  NO.     3  AT  4.23380E+01  CPS  (   2.66017E+02  RAD/SEC   ) 


FINITE  ELEMENT  MODEL  WITH  TOTAL  COMPOSITE  CASE 
12-19-91      18:06  MSC/PAL  2   V  4.0 

-AIM9  GEOMETRY  FILE 


Page 


MODE  NO.  1  AT  1.98019E+01  CPS  (  1.24419E+02  RAD/SEC  ) 
MODE  NO.  2  AT  2.14401E+01  CPS  (  1.34712E+02  RAD/SEC  ) 
MODE  NO.      3  AT  3.90748E+02  CPS  (   2.45514E+02  RAD/SEC   ) 
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TITLE  -  AIM9  GEOMETRY  FILE 

C 

C  DEFINE  NODES  FOR  MISSILE  BODY 

C 

NODE  3 

1  2.50,  0,  0  THROUGH  16  2.50,  337.5,  0 

545  2.50,  0,  102  THROUGH  560  2.50,  337.5,  102 

561  0,  0,  0 

NODE  22 

1,  16,  560,  1,  16 

NODE  22 

1,  16,  560,  15,  16 

NODE  3 

562  2.50,  0,  103  THROUGH  577  2.50,  337.5,  103 
594  1.25,  0,  112  THROUGH  609  1.25,  337.5,  112 

NODE  22 

562,  577,  609,  1,  16 

NODE  22 

562,  577,  609,  15,  16 

NODE  3 

610  0,  0,  113 

NODE  1 

17  0,0,-1  THROUGH  32  0,0,-1 

NODE  1 

145  0,  0,  -1  THROUGH  160  0,  0,  -1 

NODE  1 

529  0,  0,  -2  THROUGH  544  0,  0,  -2 

C  Missile  Centerline  Mass  Distribution  Nodes 
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,0 
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726 
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727 
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728 

o( 
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o, 

0,112 

C  DEFINE 

NODE 

POINTS 

C 

NODE  3 

611 
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C  DEFINE  NODE  POINTS  FOR  REAR  WINGS 

C 

NODE  3 

647  4.75,  45,  2 

648  4.75,  45,  6 

649  4.75,  45,  9 

650  4.75,  45,  12 

651  4.75,  45,  15 

652  4.75,  45,  18 

653  4.75,  45,  21 

654  4.75,  45,  24 

655  7.88,  45,  21 

656  11.50,  45,  18  36 

657  11.50,  45,  15 
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12 

704 
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315, 

9 
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6 

706 
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c 

C  ENTER  MATERIAL  PROPERTIES         /WcJ     T£*JbL£    Sr&£*WU 

c    njcoa^s     snort     m a4s  ast*  po&o*    i^itc 

MAT  30.0E6,  12.0E6,  7.76E-4,  0.25,  30.0E3 

C   " 

C  DEFINE  PLATE  TYPE  FOR  MISSILE  BODY 

C 

QUAD  1,  1,  0.25 

GENERATE   CONNECTS    1,    16,    560,    1,    16 

3ENERATE   CONNECTS   1,    16,    560,    15,    16 
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GENERATE 

CONNECTS  562, 

577, 

609, 

1, 

16 

GENERATE 

C 

C  DEFINE 

CONNECTS  562, 

577, 

609, 

15, 

16 

REAR 

MISSILE 

PLATE 

AND 

CONNECTIVITY 

TRI 

1, 

»  2, 

,  0.25 

CON 

1 

TO 

2  TO 

561 

CON 

2 

TO 

3  TO 

561 

CON 

3 

TO 

4  TO 

561 

CON 

4 

TO 

5  TO 

561 

CON 

5 

TO 

6  TO 

561 

CON 

6 

TO 

7  TO 

561 

CON 

7 

TO 

8  TO 

561 

CON 

8 

TO 

9  TO 

561 

CON 

9 

TO 

10  TO  561 

CON 

10  TO  11  TO  561 

CON  11  TO  12  TO  561 

CON  12  TO  13  TO  561 

CON  13  TO  14  TO  561 

CON  14  TO  15  TO  561 

CON  15  TO  16  TO  561 

CON  16  TO  1  TO  561 

C 

C  DEFINE  PLATE  CAP  CONNECTIVITY  OF  NOSE  CONE 

C 

CON  594  TO  595  TO  610 

CON  595  TO  596  TO  610 

CON  596  TO  597  TO  610 

CON  597  TO  598  TO  610 

CON  598  TO  599  TO  610 

CON  599  TO  600  TO  610 

CON  600  TO  601  TO  610 

CON  601  TO  602  TO  610 

CON  602  TO  603  TO  610 

CON  603  TO  604  TO  610 

CON  604  TO  605  TO  610 

CON  605  TO  606  TO  610 

CON  €06  TO  607  TO  610 

CON  607  TO  608  TO  610 

CON  608  TO  609  TO  610 

CON  609  TO  594  TO  610 

C 

C  CONNECT  BODY  TO  NOSECONE 

C 

QUAD  1,  lr  0.25 

CON  545  TO  562  TO  563  TO  546 

CON  546  TO  563  TO  564  TO  547 

CON  547  TO  564  TO  565  TO  548 

CON  548  TO  565  TO  566  TO  549 

CON  549  TO  566  TO  567  TO  550 

CON  550  TO  567  TO  568  TO  551 

CON  551  TO  568  TO  569  TO  552 

CON  552  TO  569  TO  570  TO  553 

CON  553  TO  570  TO  571  TO  554 

CON  554  TO  571  TO  572  TO  555 

CON  555  TO  572  TO  573  TO  556 

CON  556  TO  573  TO  574  TO  557 

CON  557  TO  574  TO  575  TO  558 

CON  558  TO  575  TO  576  TO  559 

CON  559  TO  576  TO  577  TO  560 
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CON  560  TO  577  TO  562  TO  545 

MAT  30.0E6,  12.0E6,  7.76E-5,  0.25,  30.0E3 

BEAM  3,  0.5,  0.25 

C 

C  CONNECTING  CENTERLINE  NODES  TO  EACH  OTHER  AND  TO  SURFACE  NODES 

CON  561  TO  707 

C 

CON  707  TO  708 

CON  708  TO  709 

CON  709  TO  710 

CON  710  TO  711 

CON  711  TO  712 

CON  712  TO  713 

CON  713  TO  714 

CON  714  TO  715 

CON  715  TO  716 

CON  716  TO  717 

CON  717  TO  718 

CON  718  TO  719 

CON  719  TO  720 

CON  720  TO  721 

CON  721  TO  722 

CON  722  TO  723 

CON  723  TO  724 

CON  724  TO  725 

CON  725  TO  726 

CON  726  TO  727 

CON  727  TO  728 

CON  728  TO  729 

CON  729  TO  610 

CON  707  TO  33 

CON  707  TO  34 

CON  707  TO  35 

CON  707  TO  36 

CON  707  TO  37 

CON  707  TO  38 

CON  707  TO  39 

SON  707  TO  40 

TON  707  TO  41 

20N  707  TO  42 

20N  707  TO  43 

SON  707  TO  44 

SON  707  TO  45 

ZON  707  TO  46 

:ON  707  T047 

:0N  707  TO  48 

TON  708  TO  65 

:ON  708  TO  66 

:ON  708  TO  67 

:ON  708  TO  68 

:ON  708  TO  69 

:ON  708  TO  70 

:ON  708  TO  71 

:ON  708  TO  72 

:ON  708  TO  73 

:ON  708  TO  74 

:ON  708  TO  75 

:ON  708  TO  76 
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MAT  30.0E6,  12.0E6,  7.76E-4,  0.25,  30.0E3 
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TO 

514 

CON 

726 

TO 

515 

CON 

726 

TO 

516 

CON 

726 

TO 

517 

CON 

726 

TO 

518 

CON 

726 

TO 

519 

CON 

726 

TO 

520 

CON 

726 

TO 

521 

CON 

726 

TO 

522 

CON 

726 

TO 

523 

CON 

726 

TO 

524 

CON 

726 

TO 

525 

CON 

726 

TO 

526 

CON 

726 

TO 

527 

CON 

726 

TO 

528 

CON 

727 

TO 

529 

CON 

727 

TO 

530 

CON 

727 

TO 

531 

CON 

727 

TO 

532 

CON 

727 

TO 

533 

CON 

727 

TO 

534 

CON 

727 

TO 

535 
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CON  727  TO  53  6 

CON  727  TO  537 

CON  727  TO  538 

CON  727  TO  539 

CON  727  TO  540 

CON  727  TO  541 

CON  727  TO  542 

CON  727  TO  543 

CON  727  TO  544 

CON  728  TO  545 

CON  728  TO  546 

CON  728  TO  562 

CON  728  TO  563 

CON  728  TO  564 

CON  728  TO  565 

CON  728  TO  566 

CON  728  TO  567 

CON  728  TO  568 

CON  728  TO  569 

CON  728  TO  570 

CON  728  TO  571 

CON  728  TO  572 

CON  728  TO  573 

CON  728  TO  574 

CON  728  TO  575 

CON  729  TO  594 

CON  729  TO  595 

CON  729  TO  596 

CON  729  TO  597 

CON  729  TO  598 

CON  729  TO  599 

CON  729  TO  600 

CON  729  TO  601 

CON  729  TO  602 

CON  729  TO  603 

CON  729  TO  604 

CON  729  TO  605 

CON  729  TO  606 

CON  729  TO  607 

CON  729  TO  608 

CON  729  TO  609 

C 

MAT  30.0E6,  12.0E6,  7.76E-4,  .25,  30.0E3 

C  FIRST  FRONT  FIN 

C 

QUAD  1,  1,  0.35 

CON  515  TO  531  TO  613  TO  612 

QUAD  1,  1,  0.10 

CON  611  TO  612  TO  617  TO  618 

CON  S12  TO  613  TO  616  TO  617 

CON  613  TO  614  TO  615  TO  616 

TRI  1,  2,  0.10 

CON  619  TO  618  TO  617 

CON  619  TO  617  TO  616 

CON  619  TO  616  TO  615 

C 

C  SECOND  FRONT  FIN 

C 

QUAD  1,  1,  0.35 
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CON  519  TO  535  TO  622  TO  621 
QUAD  1,  1,  0.1 

CON  620  TO  621  TO  626  TO  627 
CON  621  TO  622  TO  625  TO  626 
CON  622  TO  623  TO  624  TO  625 
TRI  1,  2,  0.1 
CON  628  TO  627  TO  626 
CON  628  TO  626  TO  625 
CON  628  TO  625  TO  624 
C 

C  THIRD  FRONT  FIN 
C 

QUAD  1,  1,  0.35 
CON  523  TO  539  TO  631  TO  630 
QUAD  1,  1,  0.1 

CON  629  TO  630  TO  635  TO  636 
CON  630  TO  631  TO  634  TO  635 
CON  631  TO  632  TO  633  TO  634 
TRI  1,  2,  0.1 
'CON  637  TO  636  TO  635 
CON  637  TO  635  TO  634 
CON  637  TO  634  TO  633 
C 

C  FOURTH  FRONT  FIN 
C 

QUAD  1,  1,  0.35 
CON  527  TO  543  TO  640  TO  639 
QUAD  1,  1,  0.1 

CON  638  TO  639  TO  644  TO  645 
CON  639  TO  640  TO  643  TO  644 
CON  640  TO  641  TO  642  TO  643 
TRI  1,  2,  0.1 
CON  646  TO  645  TO  644 
CON  646  TO  644  TO  643 
CON  646  TO  643  TO  642 
C 

C  DEFINE  REAR  WING  MATERIAL  PROPERTIES 
C 

MAT  10E6,  4E6,  2.59E-4,  0.25,  10E3 
C 

C  DEFINE  REAR  WING  PLATE  TYPE  AND  CONNECTIVITY 
C 

C  FIRST  REAR  WING 
C 

QUAD  1,  1,  0.5 
CON  19  TO  35  TO  648  TO  647 
CON  647  TO  648  TO  660  TO  661 
CON  35  TO  51  TO  649  TO  648 
CON  648  TO  649  TO  659  TO  660 
CON  51  TO  67  TO  650  TO  649 
CON  649  TO  650  TO  658  TO  659 
CON  67  TO  83  TO  651  TO  650 
CON  650  TO  651  TO  657  TO  658 
CON  83  TO  99  TO  652  TO  651 
CON  651  TO  652  TO  656  TO  657 
CON  99  TO  115  TO  653  TO  652 
CON  652  TO  653  TO  655  TO  656 
CON  115  TO  131  TO  654  TO  653 
TRI  1,  2,  0.5 
CON  654  TO  131  TO  147 
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CON  655  TO  653  TO  654 

C 

C  SECOND  REAR  WING 

C 

QUAD  1,  1,  0.5 

CON  23  TO  39  TO  663  TO  662 

CON  662  TO  663  TO  675  TO  676 

CON  39  TO  55  TO  664  TO  663 

CON  663  TO  664  TO  674  TO  675 

CON  55  TO  71  TO  665  TO  664 

CON  664  TO  665  TO  673  TO  674 

CON  71  TO  87  TO  666  TO  665 

CON  >665  TO  666  TO  672  TO  673 

CON  87  TO  103  TO  667  TO  666 

CON  666  TO  667  TO  671  TO  672 

CON  103  TO  119  TO  668  TO  667 

CON  667  TO  668  TO  670  TO  671 

CON  119  TO  135  TO  669  TO  668 

TRI  1,  2,  0.5 

CON  669  TO  135  TO  151 

CON  670  TO  668  TO  669 

C 

C  THIRD  REAR  WING 

C 

QUAD  1,  1,  0.5 

CON  27  TO  43  TO  678  TO  677 

CON  677  TO  878  TO  690  TO  691 

CON  43  TO  59  TO  679  TO  678 

CON  678  TO  679  TO  689  TO  690 

CON  59  TO  75  TO  680  TO  679 

CON  679  TO  680  TO  688  TO  689 

CON  75  TO  91  TO  681  TO  680 

CON  680  TO  681  TO  687  TO  688 

CON  91  TO  107  TO  682  TO  681 

CON  681  TO  682  TO  686  TO  687 

CON  107  TO  123  TO  683  TO  682 

CON  682  TO  683  TO  685  TO  686 

CON  123  TO  139  TO  684  TO  683 

TRI  1,  2,  0.5 

CON  <684  TO  139  TO  155 

CON  685  TO  683  TO  684 

C 

C  FOURTH  REAR  WING 

C 

QUAD  1,  1,  0.5 

CON  31  TO  47  TO  693  TO  692 

CON  692  TO  693  TO  705  TO  706 

CON  47  TO  63  TO  694  TO  693 

CON  693  TO  694  TO  704  TO  705 

CON  63  TO  79  TO  695  TO  694 

CON  694  TO  695  TO  703  TO  704 

CON  79  TO  95  TO  696  TO  695 

CON  695  TO  696  TO  702  TO  703 

CON  95  TO  111  TO  697  TO  696 

CON  696  TO  697  TO  701  TO  702 

CON  111  TO  127  TO  698  TO  697 

CON  697  TO  698  TO  700  TO  701 

CON  127  TO  143  TO  699  TO  698 

TRI  1,  2,  0.5 

CON  699  TO  143  TO  159 
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CON  700   TO    698   TO    699 
C 

C   DEFINE  ATTACHMENT   POINTS    OF  MISSILE   TO   LAUNCHER 
C 


ZERO  1 

TA  85,  213,  373 

RY  85,  213,  373 

RZ  85,  213,  373 

C 

C  DEFINE 

C 

C  MASS  C 

LUMPED  MASSES  FOR  SEPARATE  MISSILE  COMPONENTS 

IISTRIBUTION  OF  THE  MOTOR 

MASS  707, 

0.02387 

MASS  708, 

0.02387 

MASS  709, 

0.02387 

MASS  710, 

0.02387 

MASS  711, 

0.02387 

MASS  712, 

0.02387 

MASS  712, 

0.02387 

MASS  713, 

0.02387 

MASS  714, 

0.02387 

MASS  715, 

0.02387 

MASS  716, 

0.02387 

MASS  717, 

0.02387 

C  MASS  DISTRIBUTION  OF  THE  WARHEAD 

MASS  718, 

0.02133 

MASS  719, 

0.02133 

MASS  720, 

0.02133 

MASS  721, 

0.02133 

MASS  722, 

0.02133 

MASS  723, 

0.02133 

C  MASS  DISTRIBUTION  OF  THE  GUIDANCE  AND  CONTROL 

MASS  724, 

0.00583 

MASS  725, 

0.00583 

MASS  726, 

0.00583 

MASS  727, 

0.00583 

C  MASS  DISTRIBUTION  OF  THE  DETECTOR 

MASS  728, 

0.03365 

MASS  729, 

0.03365 

END 
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